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Abstract 

We consider isospin correlations of pions produced in a relativistic nuclear collision, 
using an effective theory of the chiral order parameter. Our theory has (1+1) Lorentz 
invariance as appropriate for the central rapidity region. We argue that in certain 
regions of space correlations of the chiral order parameter are described by the fixed 
point of the (1+1) WZNW model. The corresponding anomalous dimension deter- 
mines scaling of the probability to observe a correlated cluster of pions with the size 
of this cluster in rapidity. Though the maximal size of clusters for which this scaling 
is applicable is cut off by pion mass, such clusters can still include sufficiently many 
particles to make the scaling observable. 
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The Centauro events observed in cosmic ray experiments [![] stimulated interest 
in unusual states of hadronic matter 0. A now popular idea is that of disoriented 
chiral condensate P, ^, one assumes that in course of a relativistic nuclear col- 
lision, relatively large domains of nuclear matter may emit pions coherently. An 
experimental signature of such correlated domains would be that certain regions in 
the angle-rapidity plane are overpopulated by pions of the same isospin, say, there 
are regions containing mostly vr^'s without many of vr^'s, or vice versa. There is a 
proposal to look for such events at the SSC 0]. These ideas are now discussed not 
only in connection with leading-particle effects but also for the central rapidity region. 
Rajagopal and Wilczek consider formation of correlated domains at the second- 
order QCD phase transition. They argue that unless temperature drops very fast, the 
second-order phase transition is unlikely to lead to large correlated domains because 
the growing correlation length is cutoff by the pion mass. On the other hand, for a 
sudden quench from high to low temperatures, they find evidence for amplification of 
long-wavelength pion modes. 

The results of ref.0 are obtained within a (3+l)-dimensional effective theory. 
One may wonder whether this is always an adequate dimensionality. Indeed, the 
analysis of experimental data p, ^, in particular those from the CERN SPS collider 
[§], indicates that the distribution of particles in the central rapidity region possesses 
(1+1) Lorentz invariance with respect to the longitudinal coordinate z and time t. 
This means that the distribution looks alike in all reference frames obtained from the 
center-of-mass (cm.) frame by modest Lorentz boosts in the longitudinal direction. 
On the other hand, the nuclear matter produced in a collision does not maintain 
either (3+1) Lorentz invariance or the three-dimensional rotation symmetry. We will 
argue that in these circumstances, the most long-range part of isospin correlations is 
described by a (l+l)-dimensional effective theory. We also discuss a cutoff role of 
the pion mass in this case. Our analysis is simpler for the central rapidity region and 
we concentrate on this region in what follows. This does not exclude a possibility of 
similar correlation phenomena at larger rapidities. 

We need to know statistical weights of various chiral configurations in the regions 
of physical space producing the central rapidity particles. To maintain the (1+1) 
Lorentz symmetry, these statistical weights should be determined by a (1+1) Lorentz- 
invariant effective action. For low-momentum low-frequency modes, the essentially 
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unique choice is a variant of the non-hnear sigma modeL[| Since we do not have to 
maintain the full (3+1) Lorentz invariance in this effective theory, more terms become 
allowed in the chiral lagrangian. In particular, now there are two terms of the second 
order in derivatives with respect to z and t. They form the action 

/ 



where = z,t, u = exp(i t^tt"//), tt" are the three pion fields, are the Pauli 
matrices and x± are coordinates in the transverse plane. We will discuss the region of 
transverse integration below. The quantity r2 is the two-dimensional Wess-Zumino 
term |10[; it is compatible with both Lorentz and chiral symmetries. For small vr", 
the Wess-Zumino term has the form 

dzdte^"'e''^''7T^d^n^d^7T^ . (2) 



In addition to (0), there are terms containing derivatives with respect to x±; the one 
lowest in derivatives is 

- J d'^xKap Trdaud/3U'^ , (3) 

where a,P = x±i,x±2 and K^/s is some 2x2 tensor. Note that due to the (1+1) 
Lorentz invariance, there are no second-order in derivatives terms that would mix 
longitudinal and transverse indices. 

Because we want to study multiparticle, and in this sense long-range, correlations, 
we are interested in the values of the couplings /, g and Kajs attained after an appro- 
priate renormalization. A suitable renormalization procedure in our case is the scaling 
z Xz and t ^ \t leaving the transverse coordinates x±_ unchanged. Note that the 
longitudinal and transverse couplings can renormalize differently. By simple counting 
of dimensions, the transverse coupling is important when K^js ^ /x^, where /i is a 
normalization mass. At Kap <^ /i^ the scaling laws of the longitudinal theory (|l|) that 
describes the (l+l)-dimensional fluctuations renormalizing the couplings, are those 
of the (1+1) WZNW model (see below). The chiral field u scales as m oc for 7 

we take the level k = 1 value 7 = 1. The effective coupling K{^) then scales as /i"^ 
and therefore decreases in infrared. Since K decreases slower than /x^, at sufficiently 
small fi the dimensionless ratio K/ /i^ becomes large and renormalization becomes 

t Unlike ref.|^, we do not use a linear sigma model because we expect our system to be disordered 
by fluctuations of the phases of the order parameter, rather than of its magnitude. 
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complicated. In the real world, however, the accessible values of /i are bounded from 
below by the mass of pion. We will see that the renormalized value of pion mass to 
be used in our case is m ^ 73 MeV. The renormalized value of K aX fi = m is 

K^p{m)=K^p{K)[j-^ = (26 MeV)25,^ , (4) 

where we assume that the cutoff is A = 1 GeV, and Kap^Is.) = f^Sa/s, fir = 95 MeV. 
Since K{m) <^ m^, we may approximate the behavior of the system at longitudinal 
distances smaller than by the theory at the fixed point Kaf3 = 0. 

At Ka/B = fluctuations of the order parameter in different regions of the trans- 
verse plane decouple from each other and form infinitely many one-dimensional sys- 
tems. It turns out that long-range correlations in rapidity develop only in some 
regions of the transverse plane and only for a particular transverse size. We call this 
favorable size d and will estimate it shortly to be d ~ 0.8 fm. Correlations for other 
transverse sizes are short-ranged. The size d has the significance of an infrared fixed 
point: coarse graining in proper time aligns the order parameter in certain transverse 
regions over the transverse size d. Integrating eq.(|I]) over a transverse region of this 
size gives a (H-l)-dimensional theory 

S = ^l dzdt Tid^ud^u-^ + kT2 . (5) 

This is the WZNW model [|, Here A = v^/(/rf) is a dimensionless coupling 
and k is an apriori quantized, integer coefficient. Let us temporarily neglect the mass 
effects. The theory (|^) has an infrared stable fixed point at = 47r/|fc| []TU|, [TT[]. This 
gives a value of the transverse scale d for which the rapidity correlations are most 
long-ranged. 

There may be some dependence of / and d on the energy but it need not be dramatic. 
For estimates, we take / = = 95 MeV. We now have to determine the integer k. 

The actual presence of the Wess-Zumino term in the effective theory (|I]), though 
allowed by the symmetries we discussed so far, is not at all obvious. This term 
breaks G-parity of pions, so it cannot happen in a system consisting of pions alone. 
In the real system, however, there are other degrees of freedom, and G-parity can 
be broken in interactions of pions with them. For quantized coefficients like k, the 
effects of heavier particles are not necessarily suppressed by inverse powers of their 
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masses. The next lowest degree of freedom in nuclear matter are kaons. In the (3+1) 
dimensional SU{3) (that is including kaons) chiral lagrangian, G-parity of pions is 
broken by the (3+1) Wess-Zumino term [|1^. For small fields, the (3+1) Wess-Zumino 
action has the form 

TV r 

2407r2 J ti u p a , \ J 

where A*",, = 3 is the number of colors, and the 3x3 matrix 11 is a parametrization 
of an SU{3) matrix U, U = exp(zn). We may write U = a + A where a includes 
only pions and A includes only kaons and rj. Eq.([7p plays the role of a microscopic 
action with respect to our effective theory. We have to make a reduction of eq. (^ to 
the fewer fields and fewer dimensions. This means that we need a portion of eq.(|^ 
where z and t derivatives are applied only to pions but not to kaons or t]. When 
the transverse integration is restricted to a correlated region of size d, the required 
portion of r4 factorizes, 

r4 ^ -\r,Q , (8) 

where r2 is the (1+1) Wess-Zumino term (^ and 

Q--^J d^^e'^^TiPd^AdpA , (9) 

P being a projector 3x3 matrix, P = diag(l, 1, 0). The t] field does not contribute 
to eq.(|^), so the matrix A in this equation may be regarded as made purely of kaons. 
Kaons span the projective space SU{3)/ {SU{2) x f/(l)) ~ CP^, and Q is the integral 
over the correlated region of the topological density of mappings from the transverse 
plane into CP^ |l^. Eq.(|^) is its small-field limit. Therefore, the (3+1) Wess-Zumino 
action produces a non-zero (1+1) Wess-Zumino action only in the regions where the 
transverse dependence of kaon fields is topologically non-trivial. These topologically 
non-trivial configurations of kaons need not be stable; we assume that they can live 
sufficiently long to give rise to observable effects. Since the integration in eq.(|^) is 
over a finite transverse region, Q need not be integer. The minimal value of Q for 
which our dimensional reduction is consistent is Q = ±4/3, corresponding to A; = +1 
in eq.(^). Then, for the transverse scale d, eq.(|^) gives d ^ 0.8 fm. The inverse of this 
determines a typical transverse momentum of pions produced in correlated regions, 
d-^ ~ 240 MeV. 

At the critical point, that is for the favorable transverse size, the longitudinal 
correlations in the theory (P) are scale-invariant. The (unnormalized) probability 
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P{z, z'; t, t') to have the order parameter u ahgned over a distance from z to z' for a 
time from t to t' is proportional to a two-point correlation function, 

P{z,z'-t,t') oc {u{z,t)u-\z'X)) oc [{z-z'f -{t-t'f]-^'^ . (10) 

The anomalous dimension 7 is known [^], 7 = 3/(A; + 2); for A; = 1, 7 = 1. We 



are interested in the probability (|T0|) at the value of proper time r = ri at which 
the one-dimensional expansion changes into three-dimensional and the system rapidly 
hadronizes. This time is proportional to nuclear radius; for nuclei of Pb or U, Bjorken 
[^] estimates it to be ri ~ 7 fm/c. Going over to the proper time r and rapidity y via 
z = Tsinhy and t = rcoshy, setting t = t' = ti, and assuming y — y' <ti I, we get 



P(y,y')^^^. (11) 

\y - yp 



Eq. ([TTl) gives a distribution in rapidity sizes for clusters of pions with a given neutral- 
to-charged ratio, as long as the mass effects are neglected. 

In actual comparison of eq. (|ll|) with experimental data, to separate the correlation 
predicted by (11) from statistical fluctuations, one has to select high-multiplicity 
events that have clusters of pions with a far from average neutral-to-charged ratio. 
For example, one selects clusters with more than 90 or less than 5 percent of vr^'s. 
Then, one makes a histogram of occurrences of such clusters in all of the selected 
events versus their sizes in rapidity. As will be seen below, for very large clusters, the 
probability is cut off by mass effects. However, for clusters whose extent in rapidity 
is smaller than a certain Ay, the histogram should show the scaling behavior of 
probability ( pUD with the universal exponent 7. 

Now let us consider the role of pion mass and the possibility to observe correlations 
eq.([TT|) in the real world. The mass cutoff factor now appearing in eq.(0) is 

exp (-m[(^ - z'Y - (t - t')T^^) ~ exp(-mri|?/ - y'\) . (12) 

The renormalized value of pion mass is determined from the scaling equation 

ni'{m)=m\A)(^^^y' . (13) 

For A = 1 GeV, m(A) = 140 MeV, 7 = 1, we find m = 73 MeV. 

Eq. (|l^) shows that the longer is the one-dimensional evolution of the system, the 
smaller is the maximal size of clusters for which the scaling law (|ll|) can be observed. 

Ay < (mri)-i (14) 
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(assuming Ay ^ 1). This conclusion itself is not specific for our particular (1+1) 
dimensional theory; it holds whenever the (1+1) Lorentz invariance is preserved. It 
is also not unsimilar to the conclusion reached by Rajagopal and Wilczek [^] in the 
three-dimensional case if we identify the time ti as the "cooling" time of the system. 
There is, however, an important distinction between the thermalized system of ref . 
and the present case. If the system is thermalized, the number of particles produced 
by the correlation volume is essentially independent of the energy of the initial collid- 
ing particles. On the other hand, the number of particles in a fixed rapidity interval 
is known experimentally to grow with energy. Selecting high-energy events with large 



multiplicities per unit rapidity, one may be able to observe correlations (|llD even if 
the interval ( p!4|) is relatively small (see below). Increasing the atomic number of col- 
liding nuclei, on the other hand, does not necessarily enhance correlations; in fact, it 
may even suppress them. Indeed, in our theory, a correlated domain is characterized 
by the fixed transverse size d ~ 0.8 fm. According to estimates based on the hydro- 
dynamic approach |p, multiplicity per unit rapidity grows with the atomic number 
A proportionally to the total transverse area, dN/dy oc A^^^. Hence, the portion 
of dN/dy produced by the correlated transverse region is independent of A. At the 
same time, ri oc A^^^, so the number of correlated particles in the rapidity interval 
(p!4D decreases as A~^^^. In these circumstances, the best place to observe the corre- 
lations ([TTD is high-energy nucleon-nucleon collisions (SSC), rather than nuclei-nuclei 
collisions. One should keep in mind, however, that using different estimates for the 
dependence of dN/ dy on A existing in the literature may change this conclusion. 
For hadron-hadron collisions, ri ~ 1 fm/c, so the correlations extend over the rapidity 
interval Ay ~ 2.2. It seems sufficient to have 15-20 particles per unit rapidity to be 
able to observe these correlations and deduce the scaling exponent. 

To summarize, we suggest that hadronic matter created in relativistic nuclear 
collisions is best described not by an effective thermal theory but by an effective theory 
possessing (1+1) Lorentz invariance. The basis for this suggestion is the observed 
Lorentz symmetry of the distribution of particles in the central rapidity region. We 
argue that in certain regions of space correlations of the chiral order parameter are 
described by the fixed point of the WZNW model. The corresponding anomalous 
dimension determines the scaling of the probability to observe an isospin-correlated 
cluster of pions with the size of this cluster in rapidity. Though the maximal size of 
clusters for which this scaling is applicable is cut off by pion mass, such clusters can 
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still include sufficiently many particles to make the scaling observable. 
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